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ABSTRACT 


In order to scale models of devices (ships, space capsules, 
integrated circuits, etc.) and phenomena (nuclear detonations, 
space-charge widening, sloshing in fuel tanks, etc.), it is 
necessary to determine nondimensional groupings of variables. A 
computer program establishes sets of invariant or nondimensional 
sets of variables and searches for an optimum set under specified 
optimization criteria. This Fortran IV program is based on 
algorithms of integer programming. Typical running times are 
2 minutes for 10^ sets of equations. 



COMPUTER PROGRAM FOR DIMENSIONAL ANALYSIS 


By A. D. Sloan and W. W. Happ 
Electronics Research Center 


SUMMARY 

Documentation, illustrative examples, running time, and 
failure diagnostics are presented for a program useful in scaling 
models by dimensional analysis. The program, written in FORTRAN 
IV, performs exclusively integer operations and calculates a 
basis for the null space of the matrix of a transformation in 
terms of a given basis. Application to the calculation of B- 
numbers is made. The program further optimizes a matrix of 
integers under specified constraints. "Tearing" (or diakoptics) 
of large systems, defined as minimal under specified constraints, 
is performed. 

The program is available either from the authors or from the 
Project COSMIC Library of computer programs at the University of 
Georgia, at Athens. 


INTRODUCTION 

A Fortran IV program has been developed for the calculation 
of Buckingham numbers from physical laws. Buckingham numbers are 
defined as follows. Given a set of variables [ (A^) k = l,...,n], 
dimensionally expressed in terms of some fixed reference dimen- 
sions, a Buckingham number or "B" number is an n-tuple of numbers 
[ (m fc ) k = l,...,n], such that the product: 


TTK) 

k=l 


m k 


is dimensionless. For example, let the reference dimensions be 
length, mass, and time (£,m,t), and the variables be E-energy, 
m-mass, and c-speed of light. Then, these variables dimensionally 
expressed in terms of the reference dimensions by: 

E = m£ 2 /t 2 
m = m 


c = t / 1 


A "B" number for these variables is (-1,1,2) for (E) ^(m)^(c) 2 
is dimensionless. Thus, "B” numbers express formulas and may be 
used to find dimensional relationships between variables. 

In a system with a large number of variables, more than one 
"B" number may result. By finding all "B" numbers, it is possible 
to determine all the ways of decomposing the large system into 
smaller ones. For instance, if along with the variables E,m,c, 
we included t and t, in the foregoing example, we could have 
three "B" numbers (-1,1,2,00) as before; expressing E = mc^, 

(-1,1, 0 , 2 ,-2) , since the units of E are and (0 , 0 ,-l , 1 , -1) 

since the units of c are Z/t. 


SCOPE 

It is assumed that (£ k ) and [ ( 1 %) k = l,...,n] are two "B" 
numbers for the variables (A k ) and t is a real number. Scalar 
multiplication and addition among "B" numbers are defined by 


t x {l ) + (m k ) = (t l v + m.) k = 1, ,n 


Observe that: 


n 


IT (A k ) 

k=l 




-( 


n 


TT(A k ) 

\k=l 




n 


m, 


TT ( A k) 

k=l 


is dimensionless since (£ k ) and (m k ) are "B” numbers. This shows 
that the "B" numbers form a vector space, with addition and scalar 
multiplication defined as before. 

The program has the following functions . 

(1) It determines a basis for the vector space of "B" num- 
bers; that is, it determines a set of "B" numbers such that they 
are linearly independent, and such that, any other "B" number can 
be written as a finite combination of the elements of the basis. 
Such a basis is called a complete set of "B" numbers. 

(2) According to criteria previously defined (refs. 1,2,3), 
the program examines other complete sets of "B" numbers and de- 
termines a set in which: 
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(a) The minimal sum of the absolute values of the 
entries 

(b) The maximal number of zero entries 
relative to all other complete sets examined . 

PROBLEM FORMULATION 


Coding the Problem 

The input of the problem consists of: 

(1) NREFDM = the number of reference dimensions 

(2) NVARIB = the number of variables 

(3) IVAR = an (NREFDM) by (NVARIB) matrix, which acts as a 
list (indexed by the second variable) of (NREFDM) by 1 matrices. 

The J'th entry in this list, or equivalently the J'th column of 
IVAR, defines the J'th variable in terms of the reference dimen- 
sions. For example, if the J'th variable was the Boltzmann con- 
stant, which has dimensions (length) (length) (mass) / (time) / (time) 
(temperature), and if the reference dimensions were (length, mass, 
time, charge, temperature), then: 

IVAR (1 , J) = 2 
IVAR (2, J) = 1 
IVAR ( 3 , J ) = -2 
IVAR (4, J) = 0 
IVAR (5, J) = -1 

(4) M = NREFDM 

(5) N = the number of variables in the formula of interest 

(6) IFORM = a list of N numbers which indicate which variables 
are in the formula of interest. The numbers refer to the second 
subscript of IVAR. 

(7) IOPT = 1 or 2 . It indicates which optimization procedure 
will be followed. 1 means that zeros will be maximized first and 
the sum of the absolute value of entries will be minimized next, 
while 2 means that the procedure occurs in the reverse order. 
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Input Format 


The first input card contains NREFDM and NVARIB in 215 for- 
mat. The next NVARIB cards contain NREFDM numbers on each in 
1415 format and form the columns of IVAR. If Q formulas are to 
be analyzed, then there should be Q sequences of cards in the 
following form: 

IOPT in 15 format 

M and N in 215 format 

IForm (=N numbers) in 1415 format. 

Note that all input data must be in integer form. 


Example of Problem Statement (Fig. 1) 

The reference dimensions are: 

(1) Length 

(2) Mass 

(3) Time 

(4) Charge 

The variables are: 

(1) Surface area 

(2) Mass of electrons 

(3) Charge of electrons 

(4) Current density 

(5) Current 

(6) Velocity 

(7) Energy 

(8) Electric potential. 

The 12 input data cards for this problem are shown in 
Figure 1. The card at the bottom of the page is the first, while 
the one at the top of the page is the last. The resultant print- 
out is shown in Figure 2. The program listing is given in the 
appendix. 

The first card indicates that NREFDM = 4, while NVARIB = 8. 
Each of the next eight cards forms a column of IVAR. The second 
card, for example, indicated that the dimensions of variable 
number 1 (= surface area) are (length) 2 , while the ninth card 
indicates that the dimensions of variable No. 8 (= electric 
potential) are (length) 2 (mass) (time) “2 (charge) “1 . The tenth card 
indicates that IOPT = 1, and the eleventh card that M = 4, N = 8. 
The last card is IFORM. This card shows which of the variables 
are to be considered; presently all eight variables are included. 
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THERE ARE FOUR REFERENCE DIMENSIONS AND EIGHT VARIABLES. 
THE ROWS BELOW ARE THE VARIABLES. 

2 0 0 0 

0 10 0 
0 0 0 1 

- 20-1 1 
0 0-11 

10-10 
2 1-20 

2 1 - 2-1 


THE FORMULA NOW BEING CONSIDERED IS DEFINED IN TERMS OF THE 
FOLLOWING EIGHT VARIABLES. 

1 2 3 4 5 6 7 8 


THE COLUMNS OF THE MATRIX BELOW ARE THE VARIABLES IN TERMS OF 
THE REFERENCE DIMENSIONS. 

2 0 0 -2 0 1 2 2 

0 1 0 0 0 0 1 1 

0 0 0 -1 -1 -1 -2 -2 

0 0 1 1 1 0 0 -1 


THERE ARE FOUR LINEARLY INDEPENDENT B VECTORS. 

WE HAVE OPTIMIZED ZEROS FIRST, THEN THE SUM. 

THE ROWS OF THE MATRIX BELOW FORM A COMPLETE SET OF B VECTORS. 
THE SUM OF THE ABSOLUTE VALUES OF THE ENTIRES IS 19 AND THE 
NUMBER OF ZERO ENTRIES IS 19. 

-1 0 - 0-1 1 -0 -0 -0 

0 0 1 0 -0 -0 -1 1 
0 - 1-0 0 -0 -2 1 -0 

-3 0 2-2 0 2 - 0-0 


Figure 2.- Computer Print-out for data given by Figure 1. 
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COMPUTER APPROACH 


Complete Sets of "B" Numbers 

The program proceeds as follows . Suppose IVAR is a p-by-k 
matrix, so that each of the k columns of IVAR expresses a variable 
in terms of the p reference dimension. A k-tuple (ni,...,n] c ) of 
numbers is a "B" number if and only if 


0 = IVAR • 


n n 


L n k 


and this equals the matrix product of IVAR and the column vector 


n. 



If we consider IVAR as a linear transformation, then the vector 
space of "B" numbers is exactly equal to the null space of IVAR, 
and in order to find a complete set of "B" numbers, it suffices 
to find a basis for this null space. This is accomplished next. 

Let e ± = (0, ,0,1,0,...,0) be a k-tuple where the 1 occurs 

in the i-th spot, i = l,...,k. Form the matrix whose i-th row 
consists of IVAR (ei) . The matrix is k by p. Row reduce this 
matrix. As a by-product of the row reduction process one obtains 
a permutation tt of (1, . . . ,k) and number aj_g such that a^ ^ 0 and 
such that the rows of the row reduced matrix are given by 


IVAR[e^ (1) ] 
a2^IVAR[e 7r ( 3 ) ] 

a 31 IVAR[e7r (1) ] 

• • • 

a ki IVAR t e ir (1) 3 


+ a 22 IVAR[e 7T ( 2 ) ] 

+ a 32 IVAR[e Tr ( 2 ) ] + a 33 IVAR[e 7r ( 3 ) ] 
+ . . . + a^IVARte^ (k) ] 
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where the first s rows are non-zero and the last k-s rows are 
zero. Let 


W 1 = e 7T(l) 
and let ^ 

w i = S a iy e ir (y) ' 1 = 

y=l 


Since the ej/s are linearly independent and since a-j^ ^ 0 , it 
follows that the Wj _ 1 s are linearly independent. Let be the 
space spanned by wi,...,w s and let W2 be the space spanned by 
w s+ i, . . . ,Wk. IVAN (Wi©W2) = IVAR (W]_) so that W2 is the null- 
space of IVAR. Let 0 be the inverse permutation to it. Then 


f a i 0 (y) J x - 1 + 1, ,k 

expresses a basis of W 2 in terms of the original one. (Note that 
the ei could have been chosen to be any k linearly independent 
k-tuples . ) 

In this way the program determines a complete set of "B" 
numbers. The input entries of IVAR must be integers, and as a 
result, the entries of the "B" vectors forming a complete set are 
integer. 


Optimization 

Next, it is determined if there are any common factors among 
the entries of the "B" numbers, and whether there is only one 
"B" number in the complete set. If there is, there is no need to 
go through the optimizing procedure. But if there is more than 
one "B" number, then an attempt to find a more optimal complete 
set is made. The procedure is as follows. 

Suppose Vj_,... f V] c is the complete set of "B" numbers which 
we have. If 

W-^ = A + ... + A k V k where Aj_/...fAj c are integers and 
A 1 ^ 0 , then W 1 ,V 2 /-.*,V^ is also a complete set. 
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In this program, A;l,...,A}. takes on the values -2, -1,0, 1,2 and 
the different W^'s so obtained are compared and an optimal one is 
chosen. Then, Vi = V^, V 2 = Wq, V 3 = V 2 , . . . , are set. 

This process is repeated k-times. There are at present two choices 
of the optimizing procedure. The first choice is to maximize the 
number of zero entries and then minimize the sum of the absolute 
value, while the second is just the reverse. The two procedures 
do not appear to be equivalent. 

After performing the foregoing operations on one formula, 
the program goes on and reads another formula. 


Performance 

The program has been tested on an IBM 7094. The running 
time for 28 formulas, each with 1,2,3, or 4 "B" numbers in its 
complete set, was about 2 minutes. The time for a formula with 
5 "B" numbers in its complete set was slightly more than 4 min- 
utes. This time includes an assembly listing and the punching of 
a binary deck. Without the listing and deck, a Univac 1108 took 
42 seconds to do the same 28 formulas. 


Electronics Research Center 

National Aeronautics and Space Administration 
Cambridge, Massachusetts, September 1968 
129-06-02-75 
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APPENDIX 


PROGRAM LISTING 


07 

M AIN - E.FM SOURCE STATEMENT - IFN(S) - 


THIS PROGRAM PERFORMS TnE FOLLOWING FUNCTION. GIVEN A 
SET OF REFERENCE DIMENSIONS. A SET OF VARIABLES EXPRESSED 
111 TLR' ! S OF THOSE DIMENSIONS AND A FORMULA IN TERMS OF SOME" 
OF THE VARIABLES. THE PROGRAM FIRST DETERMINES A COMPLEJF 

( i . e. .maximal linearly' independent) set of b-vfctors. then" 

IT EXAMINES OTHER COMPLETE sETS OF B-VECTORS AND DETERMINES 
.HICH HAS THE MINIMAL SUM OF THE ABSOLUTE VALUES OF THE 
ENTRIES AMD WHICH HAS THE M AX IMAL NUMBER OF ZEROS IM THF 
ENTRIES OF THE b-VECTORS. 

THE' INPUT CONSISTS OF ' ' 

1. NKEFO-1 = THE NUMBER OF REFERENCE DIMENSIONS. 

Z. i.VARFTS THE NUMBER OF v API ARLES 

3. IVAR = AN (NkEFDm) BY (NVAPIR) MATRIX WHICH ACTS AS A 
LIST ( INDEXED BY THE SECOND vARIAPLE) OF (NREFD 1 -’ ) BY 1 
MATRICES. THE J’TH ENTRY IN THIS LIST » OR EQUIVALENTLY 

TliE J’Trl COLU'i'i OF IVAR DEFINES THE J’TH VARIABLE IN TERMS ' 
OF THE REFERENCE. DIMENSIONS. FOP EXA-’PLF IF THE J*TH 
VARIABLE WAS THE BOLTZMANN CONSTANT WHICH HAS PIMENSIONS 
( LENGTH ) (LENGTH) (MASS) / (TIME) (TIME) ( TEMPFR ATURF ) 

And IF THE REFERENCE DIMENSIONS WERE 
( LENGTH. MASS. TIME .CHARGE. TEi PERATURE) 

T l iEN 

IVAR < 1 . J ) —2 
I VAR (2. J) =1 
IVAR (3. J)=-2 
I VAR ( 4 . J ) =0 
IVARC3. J)=-l. 

4. |.i=f .RLFDM 

5. N= The NUMBER of variables in the formula OF INTEREST. 

6. iFoPh= A LIST OF N NUMBEi-S WHICH INDICATE WHICH 
VARIABLES ARE IN THE FORMUL/, OF INTEREST. THE NUMBFRS REFER 
TO HIE SECOliU SUBSCRIPT OF IVAR. 

t. IGPT= 1 OR 2. IT INDICATES WHICH OPTIMIZATION PROCEDURE 
'< 1 LL HE FOLLO ,','Eu . 1 MEANS THAT ZEROS WILL BE MAXIMIZED 

HIRST AND THE SUM OF THE ABSOLUTE VALUE OF THE ENTRIES WILL 
1 E NiIMlMlZEU .NEXT. WHILE 2 (-VANS THAT THE PPOCEOURE OCCURS 
IN Tilt kCVERSE ORDER. 

INPJI FORMAT 

THE FIRST I iPUT CAR,. CONTAINS NREFDM AMD NVARTR IN 
-AS F oRi-iAT . THE NEXT NVARIu CARTS CONTAIN NREFDV. NUMBERS 
ON Each IN 141 S FORMAT AND FQPM THE COLUMNS OF I V A P . 
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MAIN - EFN SOURCE ST/.TEMEMT - IFN(S) 


C IF Q FORMULAS ARE TO BE ANALYZED . THEM THERE SHOULD BE 

-'t Q SEQUENCES OF CARDS IN THE FOLLOWING' FORMT. . 

C_ _ _ 

c 

C IOPT IN 15 FORMAT. _ __ 

C " M AND N IN 215 FORMAT 

C IFORN'.(=N NUMBERS) IN 1415 FORMAT 

c 

C _ 

DIMENSION mult ( 30 ) '» N0LSAVT30 .TOT. MATS2 (30.307. T‘NMATT30> 3UT" 

DIMENSION IV ( 15. 15) . NULL ( 30 * 30 ) »IVAR(5.40) .IF0RMC15) . MATSS ( 30 » 30 ) 
r for "FORMA TT 21 5 ) 

102 FORMAT ( 1415) 

f — TO 00 TORMATTIOH THERE" ARE. T2T2'5fiTREFERE]XE' DTMEN5T0NS "ANDVr2 ."IThTTTAR I AS i 


1LES. ) 

7Q0T FORMAT T34H THE” R0W"5"“BEL0W “ARE THE VARIABLES'.! 

7002 FORMAT ( ////70H THE FORMULA NOW BEING CONSIDERED IS DEFINED IN TERM 
T S ' OF THE FOLLOW ING71 2 . TlH V A RIABLES./) 

7004 FORMAT ( /8QH THE COLUMNS OF THE MATRIX BELOW ARE THE VARIABLES IN T 

TERMS"' OF" THETEFERENCED IRDJST0NS77 

READ (5.101) NREFDM.NVARIB 

WRITE (6.7000) NREFDM.NVARIB 

WRITE (6.7001) 
dotsott j=i.n\tarib 

READ (5.102) (IVAR(I.J) .I=1.NREFDM) 

’ ’ T9 0 0" WRI TE f 6 .TR2T “ ( I V AR CT7J7 .T=T7RREFTJkT ‘ 

130 FORMAT (15) 

roo READ "(5.130) IOPT ” " " ~~ 

READ (5.101) M.N 

WRITE (6.7002FN 

READ (5.102) (IFORM( J) .U=1»N) 

WRITE T6.TU2T “( IFORMT J)TJ=T.T)7 ‘ 


DO 111 J=1.N 

TGN-TFORFTrjT 

DO 111 1=1. M 

T1TT RkATT T rjT=TVAR(T7IGNT 

WRITE (6.7004) 

' “ • - DO" 7003 1=1 »M 

7003 WRITE (6.102) ( INMAT ( I » J ) ► J=1 . N ) 

DO T ir=T.N 

DO 1 J1=1.N 

r tvti rvunw 

DO 2 12=1. N 

2 IV (I2»I2)-1 

C _ 

C FIND^ A BASI S FOR TH E NULL SPACE _ 

3 CALL NULSPA (M.N. INMA T. IV . NULL.K) __ 

_C _ THE ROWS OF N ULL ARE B_VECTORS_AND_ FORM A CO M PLETE S ET. 

C _ K IS _THE NUM BER OF B VECTO RS IN A COMPLET E SET. 

C 

C I SUMS" IS THE SUM T5F""THE _ AB"S0l-UTE " VALUE" OF THE ENTRIESTTT 
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MAIN - EFN SOURCE STATEMENT - IFN(S) 


I C THE FINAL OPTIMIZED MATRIX>AND IZERS IS THE NUMBER OF ZERO! 

C - ■ ■ . “ — - 

ISUMS=0 

izeRS=o — ' ' 

c 

■ C ARE THERE ANY NON-ZERO' B“ VECTORS 

C _ 

IF (K-l ) 200 f 620 f 801 

i 200 WRITE (6» 103) 

rCT3"T0Rf'WT C38H THE ONLT Tr~YECTOR~ IS THE ZERO VECTORTT "I 

GO_ TO 100 

C IF THERE IS ONLY ONE B VECTOR THERE IS NO NEED TO GO 

C ‘ THROUGH THE ' OPTIMIZATION POO IT NET" 

C 

"620 ~ IZtRS-0 — ~~ 

ISUMS=0 

C ------ - ■ • 

C CHECK TO SEE IF ANY OF THE NUMBERS 2>3»5»7»11 DIVIDES 

— c EVERT ENTRT~ OF “A "B “VECTOR"; 


C 

Call dI vchkTi t n t null ) 

902 DO 622 J=1»N 

TSUMS=ISUM^+IABS(NULL(lr J) ) " 

MATSZdf J)=NULL(lrJ) 

if ( nullt rrjn5z?7S2r»622 — “ 

621 IZERS=IZERS+1 

622 CONTINUE - -- - - 

GO TO 307 

C ‘ ‘ ‘ ‘ 

C SAVE ORIGINAL COMPLETE SET OF B VECTORS IN MATRIX NULSAV. 

c ™ " ' ' 


801 KK=K-1 

T5'0 "900~GTXD=1 »K 


DO 9000 JXD=1»N 

9000' ~ NUl 5 A V ’( TXB 7 JX D ) =N U L L( I XD'» JXDT - 
C 

c NOW "GO TRROaOFTOP TIM IZATT ON T - 


C 

C 

TT 

C 

C 

C 


c 

c 

c 

c 

e 


WE HAVE A COMPLETE ( I .E. » MAXIMAL LINEARLY INDEPENDENT)! 

" SET OF B VECTORSV SAY V ( 1 ) » V 12 ITT. . t V C K ) . 1 

IF W ( 1 ) =A*V ( 1 ) +B*V (2)+...+C*V(K) t WHERE A»B»...»C ARE 

INTEGERS' AND A IS NOT' ZER07THEN~WTlT V f znTTT TV ( K ) 1 S 

AGAIN A COMPLETE SET OF 8 VECTORS. IN THIS PROGRAM 
A r BV . • • t C TARE ONT THE V ALUES~-ZY- 1 .0*1.2 AND - THE 
W ( 1 ) [ S A RE OBTAINED IN THIS WAY. THEY ARE COMPARED AND 
AN OPTlM AU~ON E IS CHOSEN. THEN WE SET V( 1)=V(KT» V C2)=W( 3 
V(3)=V(2)*...»V(K)=V(K-1). THIS PROCESS IS REPEATED K 
'TIMES. 


C^ 

C 

C 

C 


C 


MAJROw INDICATES HOW MANY TIMES WE HAVE GONE THROUGH THE 
' "£BOVE~PROCESS. - - 

‘ T>0 - 301 1 f ‘ MA7JK0W=T7K ‘ 
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MAIN - EFN SOURCE STATEMENT - IFN(S) 


C ^ VALUFS 01- IMF ENTRIES IN 

C THE MOST OPTIMAL B VECTOR YET EXAMINED t AND I2ERSA IS THE 

“C NUMBER OF ZEROS. ~" ” 

C 

ISUMSA-9999 

IZERSA=-1 

" ”tup=2**kk 


c 

c NOW l»E ARE DETERMINING THE A»B» ... ,Ct WHICH WE" CALL" JZEROXr 

C MULT ( 1 ) » • • • » MULT (K-l) r RESPECTIVELY^ 

c ' - ---- ■ - - 

DO 670 I ZEROX-1 * 5 

‘ IF ( IZEROX-3) 650 » 670 1 650 
650 JZEROX= IZEROX-3 

DO 3006" INUMBl=l r IUP 
INUMB1=INUMB1-1 
DO 3006 INUMB2=1 » IUP 
INUMB2=INUMB2-1 

DO" 3006 INUMB3=l"f IUP " 

INUMB3=INUMB3-1 
DO 350 JNUMB=1»KK 
ITEST=2**( JNUMB-1) 

mulTTjnumbT=o 

IF ( AND ( I NUMB 1 t I TEST ) • NE . 0 . ) GO TO 612 

GO TO 350 " " ' " 

612 MULT ( JNUMB)=1 

1 F ( ANDTTNUMB2 » ITEST ) . NE . 0 . T'MULT ( JNUMB ) =-MULT ( JNUMB ) 

IF ( AND ( INUMB3 1 1 TEST ) .NE.O. ) MULT ( JNUMB ) =2*MULT ( JNUMB ) 


DO 351 J=lrN 


351 " 


TTTJT" 


TT" 

C 


NULL t MA JROW t . ) 15 THE NEW B VECTOR WHICH WE WILL EXAMINE. 




11 = 1-1 

352" DTT353 "'J=17N 

353 NULL (MA JROW r J ) =NULL (MAJROW r J)+(MULT( ID *NULSAV( I r J) ) 
"354 CONTINUE' 


C 

1 C — CHECK "TO" "SEE IF - ANT "OF "THE - NUMBERS 2 r 3 > 5 »7Vrr "DIVIDES" 

C EVERY ENTRY OF A B VECTOR. 

J CALL DIVC HK (MAJROWf Nr NUL L) 

C CHECK OPTIMIZING CONDITIONS. DEPENDING ON IOPT rEITHER 

~Z FIRST "DETERMINE"- THE" NUMBER OF ZERO "ENTIRES AND" THEN 5UM THE 

_C ABSOLUTE VALUES OF THE ENTRIES OR VICE VERSA. 

602 IZERO=0 

I SUM=o r 


DO 9 J5=l r N 

" ' I SUM- ISUM+ TABS ( NULL ( MA JROW r J5 JT " 
IF ( NULL (MA JROW »J5) )9»8r9 
6 I Z£RO=IZER0+I 
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non 


0 


M/UN _ - EFN_ SOURCE STATEMENT - IFN(S) 

9 CONTINUE 

GO TO (723.724) .'I OPT 

723 IF ( IZERO-IZERSA ) 3006.6001.3003 
60 01 IF (ISliM-lSUMSA) 3003 »’3C0’6» 3000 

724 IF ( I SUM- I SUMS A ) 3003*6002*3006 

60D2 TFnZERO-IZERSA)3006V3W6»3I503'~ ~ " 

3003 IZERSA=I ZERO 

. IS0MSA=ISUNf - • 

DO 3004 IK4=1 . N 

'3DU4' MATSZTMAJRW. IK4T =NUULTMATJR0W . IK4? 


■" MATSZ* I5"THE”M0Sr OPTTMAU' NEW B ■ VECTOR ' TFT' ' EXAM I NED". 


3006‘"C<JnTIW£ 

670 CONTINUE 

D^BIT'JSITN 

760 NULSAVd. J)=NULSAV(K. J) 
DOTTGIT" r-2~»K" 

DO 761 J=1.N 

— • 761 TWTSST1VJ)=1TOD5AV (T.'TJ 

DO 762 1=2. KK 

D 0 “762~J=rr'N 

12=1+1 

762 NULSAV (I2»J)=MATSS(T»J) 

DO 763 J=1»N 

T6T'KIDC5A'VT27 JrSMATSZ ( MAJR0W . J ) 

ISUMS=ISUMS+ISUMSA 

T7eR5=TZER5+TZERSA 

3014 CONTIJNUE_ 

C 

C WE”'ARE DONE TESTIW. 


C 


U K 


104 FORMAT ( 10H THERE ARE.I3.32H LINEARLY INDEPENDENT R VECTORS.) 

GO TO ( T31T132T7T0PT 

131 WRITE (6.133) 

'133"F0RMATnP+H ' WE "HAVE - OFT I M IZ E D' - ZEROSTTRST .THEN" THE 5UMD 

GO TO 4006 

132" WRITE (6.134) - -- •- - — 

134 FORMAT ( 44H WE HAVE OPTIMIZED THE SUM FIRST. THEN ZEROS.) 

4GU6" "WRITE T57T07) ~ — “ ' ' ‘ 

107 F0KMAT(63H -THE ROWS OF THE MATRIX BELOW FORM A COMPLETE SET OF B V 

iectorst) 

WRITE (6.120) ISUMS. IZERS 

120 'F0RMATT49FT THE - SUM OF THE ABSOLUTE" VALUES OF THE EN T R I ES 1S.I3.34H 
1 AND THE NUMBER OF ZERO ENTRIES IS. 13) 

CO 3T5~T9=i.K‘ - - . 

315 WRITE (6.106) (MATSZ( 19. J9) . J9=l .N) 

T06 FORMA TTT5T5) ~~ ~ 

GO TO 100 

End ' ' ■ 
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SI 


EFN SOURCE STATEMENT - IFN(S) 


~TT 


SUBROUTINE BABE CM » N » TNMAT » MATOCTT »K ) 

DIMENSION INMAT (30*30) .MATOUT (30 » 30 > » IC0NS2(15) » MATSAV ( 30 * 30 ) 

C ‘ ~50BT?0UTTNEr~BA5Er TTOS THEFOLUO'dNG FUNCTION. GIVEN 7T 'VECTOR 

C SPACE V WITH BASIS ( VI .V2. • • . .VN) AND A SET OF VECTORS IN V (W1»W2» 

C •' V.'7»WK) » BASE "DETERM IHES-T7K=THE"Tn MENS ION OF 'W.WHERE' TJ - ISTHE" - 

C SUBSPACE GENERATED BY THE WI’S.AND 

C" “• -2 ."K - VECTORS - TN "¥ .TTTTERMS ' OF THE~ORTOTN'AE 

C WI*S. WHICH ARE A BASIS FOR W. 

C 'THE ~rRPUT _ CONSTSTS‘'OF^~I.N=THE DIMENSION OF \T 

C 2.M=TH£ NUMBER OF VECTORS IN THE GIVEN 

C BET ► AND 

c 3.INMAT=THE M by N MATRIX WHICH 

C ’ EXPRESSES" THE' WT* 5 TN TERMS 'OF TRT'VT’5. 

C THE OUTPU T CO NS ISTS OF- I»K=THEDIMENS^ ° F ..^* A NP 

C FIRST K ROWS DETERMINES THE BASIS OF W IN TERMS OF THE VECTORS IN 

Z~ THE' GIVEN BETT 

C NOTE-THE INPUT MATRIX MUST HAVE INT EGRAL ENTRIES. 


DO 61 I = 1.M 
DO 61 J=17M 


MATSAV(I»J)=INMATd.J> 
IF ( I-J) 1.60.1 
1 MATOUTd . J)=0 
6O~T0 ST 


60 MATOUT (I »J>=1 


C 

£ 

C 

C 

z 

c 


IROW=l 



PROCESS. IT ALSO INDICATES HOW MANY VECTORS IN THE BASIS 


HAVE" 


FIND A NON-ZERO VECTOR 


DO 50 J=1.N " 

DO 50 1=1. M 

— — tfctnmxttt rj) tsttsotst 

51 JSTART=J 

rrrr-MT30T33.33 
30 ISTART=I 

SO TO 31 

_ 50 CONTINUE 

RETURN 

JTdPLU5=TBTffRT+I 


DO 32 1= I PLUS . M 
DO ~32 J=TTn 


IFdNMATd.J) >52.-32.52 

'3'2'COnTINJE ^ 


33 K=1 

RETURN ' ~ 

52 DO 15 J=JSTART.N 
DO 3 7=IR0W»M 
IF (INMATd.J) >2.3.2 

~2"r5woF=I — 
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IFN(S) 


SI _ - EFN SOURCE STATEMENT - 

GO TO 4 
5 CONTINUE 
GO TO 15 


INTERCHANGE ROWS OF INMAT. 


4 DO 5 J1=1*N 
‘SAVE=INMflTCTROW» J1 ) 
INMAT(IROW*Jl>=INMAT(ISWOP*Jl) 

5 INMAT( ISWOP* J1>=SAVE 

INTERCHANGE CORRESPONDING ROWS OF MATOUT. 
DO 6 11=1. M 
SAVE=MATOUT ( IROW .11) 

MATOUT ( I ROW *11) =MATOUT ( ISWOP *11) 

6 MATOUT C ISWOP* II )=SAVE 
9 IR0W1=IR0W+1 

IC6NSi=INMAT(IROW*J) 

ROW REDUCE INMAT. 

CHANGE MATOUT ACCORDINGLY. 


DO 66 I3=IR0W1*M 
IC0NS2(I3)=INMAT(I3*J) 

IF(IC0NS2(I3) >65*66*65 

65 DO 10 J3=1*N ' “ 

10 INMAT(I3* J3) = (IC0NS1*INMAT(I3*J3> >-(IC0NS2(I3)*INMAT(IR0W* J3> ) 

66 CONTI NUET 

DO 68 I4=IR0W1*M 

IF(IC0NS2(I4> >67*68*67 

67 DO 13 J4=1*M 

TTMATOUT (14 ► J4 ) = ( IC0N51*MAT0UT C 14 * J4) J - ( IC0N52 (14) *MATOUT ( IROW * J4) > 

68 CONTINUE 

IF(J-N>40*I9VW 

40 IF ( I R0W1~M) 41 *16*19 

41 JONE=J+l - - 

DO 42 JTRY=JONE*N 

'UCT _ 42' ITRT=rR 0 Wl * M 

IF l INMAT ( ITRY * JTRY )) 14*42* 14 

42 CONTINUE 

GO TO 19 

14 IR0W=IROW+l " 

15 CONTINUE 

go 10 iy - — 

16 DO 17 15=1 *N 

IF (INMAT ( lR0¥l*TSrn 18* 1T»TIB " 

17 CONTINUE 

19 K=IROW ^ 

GO TO 71 

— IB “K=TROWT — 

71 DO 70 I6=1*M 

. D O *70 I7=17N 

70 INMAT (16* I7)=MATSAV (16*17) 

i 'RETURN 

END 
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S2 


EFN SOURCE STATEMENT 


IFN(S) 


Ql/, 

i 


c 

c 

c 

■ c 
c 

~T~ 

c 

~ZT 

c 

c 

c 

c 

c 

c 

c 

c 

~zr 

c 

c 

c 

“C - 

c 

c 

c 

~cr~ 

c 

c 

“C - 

c 

“C - 

c 

"C — 

c 

“X 

c 

c 

c 

c 

X — 

c 

X” 

c 

c 


SUBROUT I NE NULSPA ( M » N » L » I V » NULL • K ) 

WE ARE GIVEN TWO VECTOR SPACES V AND W WITH BASIS V(l)»... 

V (N) AND W(l) t . . . »WCM) .RESPECTIVELY. AND A LINEAR TRANSFORMA- 

TION L FROM V INTO W. SUBROUTINE NULSPA DETERMINES A BASIS FOR 

THE NULL SPACE OF L IN TERMS OF A GIVEN BASIS V. 

VTTT.....V(N> IS THE USUAL BASIS. NAMELY VTTTrTTT. . . . » 0 » 

1.0.... .0) WHERE THE 1 OCCURS IN THE I*TH LOCATION. 


THE INPUT CONSISTS OF- 


1. TWO INTEGERS M AND N, 


_2.AN_N BY N MATRIX IV WHOSE COLUMNS AR E A BASI S OF V. AND 
3. AN M BY N MATRIX WHICH CONRRESPONDS TO L IN THE GIVEN BASIS. 


TfflTWTPUT CONSlSTS”OF- 

_ r.K =' THE DIMENSION OF THE NULL SPACE OF L . AND 

~2.NULL= A "K ST N“MATRIX WHOSE U*TH ROW <Jl». ..UN) REPRESENTS " 
THE J’TH VECTOR = U1*V<1) +...+JN*V(N) OF THE BASIS OF THE j 
TTULL SPACE. " 


THE METHOD. 

LET V CD BE THE I»TH COLUMN VECTOR OF IV. STEP 1 
CON SISTS OF FOR MING THE MATRIX IRANGE WHOSE I’TH ROW CONSIS 
T5 OF HVII) )'« STEP 2' CONSISTS OF CALLING SUBROUTINE BASE 
WHICH ROW REDUCES IRANGE. IT ALSO DETERMINES THE ROWS OF 
THE ROW REDUCED MATRIX INTERMS OF THE INITIAL ROWS OF THE 
MATRIX AND EXPRESSES THE RESULT IN THE MATRIX IBSRAN. 

IDMRAN IS THE DIMENSION OF THE RANGE OF L. THE 
FIRST IDMRAN ROWS OF IBSRAN FORM A BASIS OF THE RANGE OF L 


~RS~FOLLOWS 

* WE HAVE V(l)»...rV(N> AS A BASIS OF V AND WE FORMED 
LfV(I) ) FOR I-l» ... *N« SUPPOSE IBSRAN=<A ( I#U) J . THEN THE 
I*TH ELEMENT OF A BASIS FOR THE RANGE OF L IS GIVEN BY 
AiI»l)*L(VCl))+. ;.+A(I*N)*L(ViNJ) . SINCE L IS LINEAR IT 
FOLLOWS THAT THE VECTORS A( I t 1 > *L( V( 1 )+. . .+A( I »N)*V<N> t 


for r=i> .“. ; idmran>fokm a basis for a subsface of v which is 
ISOMORPHIC » UNDER THE RESTRICTION OF L TO THIS SUBSPACEt 
TCTHt RANGE OF L. THE LAST N- IDMRAN ROWS OF THE ROW 
REDUCED IRANGE MATRIX ARE ALL ZERO. SO A ( I » 1 ) *L( V( 1 ) ) +. , .+ 
A(r»NJ*L(V(N> JSOrFOR I=IDMRAN+1 t . . . »N. THE LINEARITY OF L 
SHOWS THAT THE VECTORS A ( I » 1 ) *V ( 1 ) +. . .+A ( I »N> *V (N) r FOR 1= 
TUMR AR+T N » ARE IN THE' NULL SPACE OF L, THE ROW 
REDUCTION PROCESS SHOWS THAT THESE VECTORS ARE LINEARLY 
INDEPENDENT • HENCE THE LAST N-IDMRAN ROWS OF IBSRAN ARE A 
BASIS FOR THE NULL SPACE OF L. 


i 


DIMENSION L(30»30) » NULL (30.30) 
DTWERSTUKT I V ( 157X5 ) 
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S2_ -_EFN SOURCE STATEMENT IFN(S) 

i DIMENSION I RANGE ( 30 » 30 ) t IBSRAN ( 30 » 30 ) 

r "Do i rr=i*M 

DO 1 J1=1»N 

: r trance rji7in=o" -*• — • 


STEF 1 

DO 7 12=1 »M 

— D0"T 'J2='1TN — 

DO 7 J3=1»N 

7nrRANGE"nJ57T2T=IR7fNGET'J27T?7 + '(L ( 127 J3) *IV ( J3VJ271 

STE p--2 

CALL BASE ( N * M * I RANGE ► IBSRAN ► IDMRAN ) 

- K2N- IDMRAN 

DO 6 12=1 » K 

I3=TDnTRAN+I2 

DO 6 J2=l t N 

6 NULL Cl2 r J2T=TBSRAN (TT»U21 
RETURN 
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n r> 


u 


? 


S3 - EFN SOURCE STATEMENT - IFN(S) - 


SUBROUTINE DIVCHK (L2»N»NULL) ' ‘ ' ' 

DIMENSION NULL 1 30 » 30 ) _ _ 

DIVCHK CHECKS TO SEE IF THE NUMBERS 2»3»5»7»11 DIVIDE EVERY ENTRY 

T OF “ROW L2 IN TRE‘‘MATRIX“NULL~WHrCH HAS "N‘ ' COLUMNS’* IF - SO THEN 

C THE DIVISION IS DONE. _ 

ICNT=1 

732~ GO TO (726V727 .728 » 72977317) "TI'CNT 

' 726 L3-2 

r BO'TO 731' 

! 727 L3=3 

1 '60‘ T0“ 731' " - ■ — - - 

728 L3=5 

*“ GO TO "731 

729 L3=7 

— “GO TO 731 

730 L3=ll 

‘ "73l~0O 72'0“L4=T»N" "" - - 

IChECK = MOD (NULL (L2 f L4 ) rL3) 

inTCHEOKT72^77207 7 22 “ 

720 CONTINUE 

T50 721 L5=1»N 

721 NULL (L2 * L5) = NULL ( L2 * L5) /L3 

~GCT10'~7ZZ 

722 IF(ICNT-5)733.802>802 

733" TCnT= I CFiT +TT “ ~ 

GO TO 732 

8C2 — RETURn ■ 

END 
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